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ABSTRACT 

We analyze the Chandra X-ray spectrum obtained with the HETGS grat- 
ing spectrometer of IRAS 13349+2438, which has one of the richest absorption 
spectra of a quasar outflow. Absorption from almost all charge states of Fe 
is detected. This allows for a detailed reconstruction of the absorption mea- 
sure distribution (AMD), which we define as the continuous distribution of col- 
umn density as a function of ionization parameter. We find a double peaked 
AMD for IRAS 13349+2438 with a total (ionized) column density of N H = 
(1.2+0.3) xlO 22 cm -2 assuming solar iron abundance. For comparison, we per- 
form a similar analysis on the well studied HETGS spectrum of NGC 3783. Both 
sources feature a deep minimum in column density consistent with no absorption 
from gas at temperatures of 4.5 < logT < 5 (K). We interpret the minima as 
observational evidence for thermal instability in this temperature regime. 

Subject headings: galaxies: active - - galaxies: individual (IRAS 13349+2438, 
NGC 3783) — techniques: spectroscopic — X-rays: galaxies — line: formation 



1. INTRODUCTION 

The X-ray spectra of many active galactic nuclei (AGNs) viewed directly toward the 
central source (e.g., Seyfert 1 galaxies) show the continuum flux absorbed by numerous ab- 
sorption lines produced by ionized gas. The lines are generally shifted to shorter wavelengths, 
indicating an outflowing wind. It is not clear whether this wind can also be associated with 
the line-of-sight material responsible for the optical reddening of these sources. To that end, 
the outflow would have to entail neutral- gas or dust grains along with the highly ionized 
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wind. It is still controversial whether the winds are significant for the AGN central engine 
or for the host galaxy in terms of mass, energy, and momentum. It is possible that extreme 
outflows may even play a central role in cosmological feedback and in the metal enrichment 
of the intergalactic medium (IGM). In order to address these open questions, it is important 
to develop reliable techniques to probe the winds and to perform accurate measurements of 
their physical properties. 

The X-ray band is advantageous for AGN outflow studies not only because most of the 
outflow appears to be highly ionized, but also since the X-ray band comprises detectable 
absorption lines from the full range of charge states: From neutral up to hydrogen- like. 
Moreover, in most cases an X-ray spectrum features several lines from a single ion, which 
helps to deduce ionic column densities from saturated absorption lines. The detection of 
many charge states of a single element is key to plasma diagnostics as it allows to model the 
distribution of ionization and to measure the total column density independent of assump- 
tions regarding the relative element abundances. For this purpose, Fe ions are particularly 
usefu l as they form over seve ral orders of magnitude of the ionization parameter, defined 



(e.g., 



Kallman fc Kroliklll995h by 



L 



where L is the ionizing luminosity between 1 and 1000 Rydberg, n H denotes the hydrogen 
number density and r is the distance of the absorber from the ionizing source. Once the 
distribution of ionization is known for one element, the relative abundances of other ele- 
ments can be measured by comparison. For the absolute abundances of metals with respect 
to hydrogen, one still needs to utilize UV spectra, as the energies of atomic transitions in 
hydrogen are too low to be observed in the X-rays. Highly ionized K- and L- shell ions (1-2, 
and 3-10 electrons, respectively) of Fe appear regularly in astrophysical X-ray sources, either 
in emission or in absorption. The lower charge state s of Fe (H-26 elec trons) feature strong 
inner-shell absorption lines in the soft X-ray band (IBehar et al.ll200ll ). Each ion can pro- 
duce dozens of overlapping lines forming unresolved transition arrays (UTAs) . In terestingly, 



IRAS 13349+2438 was the first source in which Fe M-shell UTAs were observed ( iSako et al. 
200lh . 



The majority of works on the X-ray spectra of AGN outflows fit the data with a grad- 
ually increasing number of iso-£ absorption components, until the fit is satisfactory (e.g., 
Kaspi et al.ll200ll ; iKaastra et al.l 120021 ; iNetzer et al.l 120031 ) . For high-quality spectra two or 
three ionization components produce a good fit. In this work, we employ a more funda- 
mental method, which tries to reconstruct the actual distribution of the column density in 
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the plasma as a continuous function of £. We call this distribution the absorption mea- 
sure distribution ( AMD ). The present method i s sim ilar in principle to that tested by 
Steenbrugge et al.l (120051 ) and by ICostantini et al.l (120071 ). who assumed a power- law distri- 
bution. However, here we fit for the AMD directly from the measured ionic column densities 
making no a-priori assumptions regarding its shape. 

IRAS 13349+ 2438 is a bright type-1 quasar (L > 10 46 erg s" 1 ) at a redsh ift of z=0. 10764 
( jKim et al.l Il995l ) . Optical and infrared observations by IWills et al.l (119921 ) h ave suggested 



absor ption of the central source by dust. X-ray observations with ROSAT (IBrandt et al. 



19961 ) and ASCA (IBrandt et al.l 1 19971 ) did not find a significant column density of cold mate- 
rial, but sug gested the pre s ence o f ionized gas along the line of sight, perhaps with embedded 
dust grains. ISiebert et al.l (119991 ) restated the difficulty to explain the optical and X-ray ob- 
servations with the same absorber. The line-resolved RGS (Reflection Grating Sp e ctrom eter 
on board XMM-Newton) spectrum of IRAS 13349+2438 measured by ISako et al.l (120011 1 re- 
vealed broad (FWHM ~ 1400+ 300 km s _1 ) absorption lines of many ions. Two distinct 
ionization components were determined at log£ ~ and at log£ ~ 2 — 2.5 where £ is in 
c.g.s. units (erg cm s _1 ). The low- and high- £ components were found to have column 
densities of a few 10 21 cm -2 and a few 10 22 cm -2 , respectively. The good a greement between 
the low-£ column de nsity and that deduced from the optical reddening (IWills et al.lll992l ) 
led lSako et al.l (120011 ) to conjecture that the low-£ absorber might be responsible also for the 
optical reddening. In this paper, we revisit the ionized absorber of IRAS 13349+2438 using a 
long exposure with the High Energy Transmission Grating (HETGS) spectrometer on board 
Chandra performed 3 years and 9 months after the XMM-Newton observation. The spectral 
resolution of HETGS is much better than that of the RGS and the long exposure compensates 
for its inferior effective 31*6 3. 



2. DATA REDUCTION 

IRAS 13349+2438 was observed by Chandra/BETGS on 22-26 February, 2004 for a total 
exposure time of 300 ks. All observations were reduced from the Chandra archive using the 
standard pipeline software (CIAO version 3.2.1). The total number of counts in the first (plus 
and minus) orders between 2 and 20 A is 22836 for MEG (medium energy grating) and 9542 
for HEG (high energy). No background subtraction was required as the back ground level 



was n egligible. The post-pipeline procedure is similar to the method reported in lKaspi et al. 



( 120011 ). Flux spectra were obtained by dividing the count spectra by the effective area 
curves produced by the standard CIAO tool. Positive and negative first-order data of both 
the MEG and HEG instruments were combined. The spectrum was corrected for neutral 
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galactic absorption (Nh = 1.1 xlO 20 cm" 2 , iMurphy et al.lll996l ). The data were rebinned to 
36 mA bins to improve the statistics of the spectra with minimal compromise of the spectral 
resolution, as the MEG instrumental line spread function is approximately 23 mA broad 
(FWHM). 



3. SPECTRAL MODEL 



Variations of approximately 25% on time scales of 20 ks were observed during the 300 ks 
exposure of IRAS 13349+2438. Due to this low level of variability and since we are interested 
in the over all properties of the ionized absorber, henceforth, we use only the combined 
MEG and HEG full 300 ks spectrum. The present fitting procedure follows our ion-by-ion 
fitting method (IBehar et al.ll200ll . 120031 ; iHolczer et al.ll2005l ). First, we fit for the broad-band 
continuum. Subsequently, we fit the absorption features using template ionic spectra that 
include all of the absorption lines and photoelectric edges of each ion, but vary with the 
broadening (so-called turbulent) velocity and the ionic column density. 



3.1. Continuum Parameters 



The continuum spectrum of IRAS 13349+2438 can be characterized by a high-energy 
power-law and a blackbody that rises at lower energies. The present best-fit power-law 
parameters are a photon index of T e = 1.9 with a total flux in the 5-35 A range of 2.8x 10~ 3 
photons cm " 2 s~ l , to be com pared to T e = 2.2 and a flux of 2.4 xlO -3 photons cm -2 s _1 
obtained by ISako et al.l (120011 ) . The current blackbody has a temperature of kT = 105 eV 
and a 5-35 A flux of 7.8 x!0~ 3 photons cm" 2 s~ \ compared to kT = 100 eV and 4.0xl0~ 3 
photo ns cm -2 s _1 obtained by ISako et al.l (120011 ). The continuum applied by iBrandt et al. 



(119971 ) to the ASCA data produced a flux in the 2-10 keV range of 5.7x10 12 erg cm 2 

The continuum applied 



s 1 compared to the p resent value of 3.6x10 



-19 —2 —1 

erg cm s . 



by IBrandt et al.l ([1996) to the ROSAT data produced a flux in the 0.1-2.5 keV range of 



2 s" 1 . 



In short, 



1.2x10 11 erg cm 2 s 1 compared to the present value of 1.7x10 11 erg cm" 
the continuum measured here with Chandra is within ~50% of all previous o bservations. 



In pa rticular, the continuum model parameters are fairly similar to those of ISako et al. 



(120011 ). In Fig. [I] we present a direct comparison of the two grating spectra from these two 
observations showing their over all similarity. It can be seen that the hard X-ray flux is 
almost identical, while there is some increase in the soft flux. This change is reflected in the 
small difference between the two models in the blackbody normalization parameters. One 
might notice that the depth of the O K-shell edges is also different between observations. 
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This is discussed in more detail below in §4.11 



3.2. The Ionized Absorber 

The spectrum lijiy) around an atomic absorption line % — > j can be represented by: 



Iij(v) = I (v) exp [-N ion aij(v)} 



(2) 



where Iq(v) represents the continuum intensity, <J%j{v) denotes the line absorption cross 
section for photoexcitation (in cm 2 ) from ground level % to level j. If all ions are essentially 
in the ground level. N ion is the total ionic column density towards the source in cm -2 . The 
photoexcitation cross section is given by: 



7re 
m P c 



(3) 



fij denotes the absorption oscillator strength and <f>(u) represents the Voigt profile due 
to the convolution of natural (Lorentzian) and Doppler (Gaussian) line broadening. The 
Doppler broadening consists of thermal and turbulent motion, but in AGN outflows, the 
turbulent broadening appears to totally dominate the temperature broadening. Transition 
wavelengths, natural widths and oscillator strengths were calculated using th e Hebrew Uni- 
versity Lawrence Livermore Atomic Code (HULLAC. iBar-Shalom et al.ll2001l). Mor e recen t 



and improved atomic data for the Fe M-shell ions were incorporated from 



Gu et all (120061 ) . 



The Doppler broadening can be typically characterized by a turbulent velocity ikurb, 
although for very high resolution UV spectra of AGN outflows, this approximation has been 
shown to be inadequate and velocity-dependent partial covering effects need to be taken 
into account (lArav et al.lll999l ). Due to the limited resolution, the present HETGS spectrum 
of IRAS 13349+2438 does not warrant a more sophisticated analysis than the v tur b approx- 
imation. Indeed, a turbulent velocity of Vt ur b = 640 km s _1 , corresponding to a FWHM of 
~ 1500 km s -1 , provides a good fit to the strongest absorption lines in the spectrum as 
demonstrated in Fig. [2]for the Fe +16 absorption line at 15.01 A. This specific absorption line 
is shown because it has a high oscillator strength of 2.3 and i t is relative l y unb lended. The 
value of v tur b = 640 km s -1 is consistent with that quoted by lSako et al.l (120011 ). 



Since the absorbing gas is outflowing, the absorption lines are slightly blue-shifted with 
respect to the AGN rest frame. The curren t best- fit outflow velocity is —300 ± 50 km s^ 1 . 
Errors correspond to A% 2 = 1. ISako et al.l (120011 ). on the other hand, found from the RGS 
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spectrum two velocity components, high- and low- ionization components with velocities of 
+20 ^200 km s _1 , and -420 1]Jq km s _1 , respectively. We find that a single outflow velocity 
of —300 ± 50 km s _1 actually provides a better fit to the HETGS spectrum. In Fig. [21 
the centroid of the Fe +16 line ca n be seen to be off set from the zero-velocity position by at 
least six sigma (er = 50 km s -1 ). ISako et al.l (120011 ) ascribed this line to the high-ionization 
component with v out = +20 l^oo km s_1 - The —300 km s _1 blue-shift (obtained by fitting all 
absorption lines in the spectrum simultaneously) clearly provides a better fit to the profile 
of this specific line. We ascribe the potential discrepancy with ISako et al.l (120011 ) to the 
superior spectra l resolution of the present HETGS measurement and to the large errors on 
the velocities of ISako et al.l (120011 ) . If anything, the data in Fig. [2] hint at an even higher 
velocity for this specific line, which might suggest that there is some scatter in velocity 
between the different ions. However, this effect is too weak for us to measure even with 
HETGS and henceforth, we assume a uniform outflow velocity of —300 ± 50 km s _1 . 

Our model includes all of the important lines of all ion species that can absorb in the 
waveband observed by HETGS. In the present spectrum of IRAS 13349+2438, we find evidence 
for the following ions: N+ 6 , all charge states of O, Ne+ 6 - Ne+ 9 , Mg+ 4 - Mg +n , Na +1 °, Si+ 5 - 
Si +13 , 21 ions of Fe, and Ni +18 . We also include the K-shell photoelectric edges for all these 
ions although their effect here is largely negligible. When fitting the data, each ionic column 
density is treated as a free parameter. A preliminary fit is obtained using a custom made 
Monte-Carlo fitting method applied to the entire spectrum. Subsequently, the final fit is 
obtained for individual ionic column densities in a more controlled manner, which ensures 
that the fit of the leading lines is not compromised. The best fit model for the present 
spectrum of IRAS 13349+2438 is shown in Fig. [3l It can be seen that the Si, Mg, Ne, and 
Fe L-shell lines are reproduced fairly well by the model. However, in the Fe M-shell region 
between 15.0 - 17.5 A, the spectrum is rather noisy and it is difficult to identify individual 
ions. For these ions, only rough estimates of the column densities can be obtained. Note 
that some lines could be saturated, e.g., the leading lines of +6 and +7 . In these cases, 
the higher order lines with lower oscillator strengths are crucial for obtaining reliable A^ on 
values. 



4. IONIC COLUMN DENSITIES 
4.1. Results 

The best-fit ionic column densities are listed in Table [TJ For the most part, the column 
densities of the Fe, Si, and Mg ions are of the order of 10 16 -10 17 cm" 2 . The column densities 
of N +6 and the Ne ions exceed 10 17 cm -2 , while those of the O ions reach ~ 10 18 cm -2 . 
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Comparing these results with those of ISako et al.l (120011 ). we find that the Fe M-shell column 
densities are consistent. The present column densities for the O K-shell io ns are much 
higher, while those of Fe L-shell are much lower than those of ISako et al.l (120011 ) . In order to 
understand this discrepancy, we went back to the RGS spectrum and compared both data 
and models with the present spectrum. We conclude that the O K-shell column densities 
indeed did vary between observations. This can be seen even in our Fig. [1] where features of 
+6 and +7 , marked by labels, are clearly deeper in the present spectrum, even though the 
soft excess continuum is higher. Conversely, for all non-oxygen ions, we find that the present 
ionic column densities obtained from fitting the HETGS spectrum actually reproduce the 
RGS spectrum fairly w ell. Other than O, we have no good explanation for the discrepancies 
with ISako etaD (bpOlh . 



Interestingly, we find a relatively high column density of neutral Fe which may not be 
part of the outflow, as manifested by the La absorption U TA at 17.53 A (see Fig. [3]). If we 
use the solar Fe/H abundance ratio (lAsplund et al.ll2005l ). we obtain a neutral gas column 
density of Nh ~ 3.5 x 10 21 cm -2 , much higher than the galactic column Nh ~ 10 20 cm -2 . 
This may indicate the presence of neutral gas or dust in the host galaxy reminiscent of the 
Wills et all fjl992h finding. 



4.2. Errors Determination 

The column density errors quoted in Table [1] were calculated by keeping an individual 
ionic column density fixed while performing a minimum - \ 2 fit with all other parameters, 
and then requiring that A% 2 = 1 for variations of that individual ionic column density, 
while all other parameters remain fixed. The more statistically rigorous method would be to 
require A% 2 = 1 following a refit of the data. However, since we supplement \ 2 fitting with a 
more controlled manual fit to the strongest lines, leading to a best fit that is not strictly the 
X 2 minimum, this method which involves automatic x 2 minimization is inapplicable here. 
Nevertheless, our simpler method gives reliable upper limits to the errors. 

For the cases in which the relative error was less than 10%, we assume an error of 10%. 
The reason is that even for the most prominent lines, we do not trust the column density 
to better than 10%. Conversely, when the relative error was more than 200%, we quote a 
nominal error of 200%. In these cases we believe the errors are not statistical. 

In §5.11 we need to use symmetrical errors. In order to get a reliable symmetrical error 
estimate, we introduced one more adjustment. We define a geometrically averaged error : 
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GTT symmetric \/^"^up ' GVT 'down- (4) 

We use geometrically averaged rather than arithmetically averaged errors, since the geomet- 
rical average gives an average in order of magnitude between the two errors, which is the 
more appropriate value to be applied as a symmetrical error. For the full AMD analysis 
(presented in §5.11) . we use the true asymmetrical errors in the manner explained below. 



5. ABSORPTION MEASURE DISTRIBUTION 
5.1. Method 

The large range of ionization states present in the absorber strongly suggests that the 
absorption arises from gas that is distributed over a wide range of ionization parameter. The 
total hydrogen column density Nh along the line of sight can therefore be expressed as an 
integral over its distribution in log£. We call this continuous distribution the Absorption 
Measure Distribution {AMD by analogy to the Emission Measure Distribution in emission- 
line spectra): 

AMD = dN H /d(\og£) (5) 

and 

Nh = J AMD d(logO (6) 
The relation between the ionic column densities N ion and the AMD is then expressed as: 

N ion = A Z [ J^ /ion (log£)d(log£) (7) 

where N ion is the measured ion column density, A z is the element abundance with respect to 
hydrogen assumed to be constant throughout the absorber, and /ion (log £) is the fractional 
ion abundance with respect to the total abundance of its element. Here, we aim at recovering 
the AMD for IRAS 13349+2438. 

As an initial approximation to be relaxed later, let us assume that each ion forms exclu- 
sively at the ionization parameter ^ ma x where its fractional abundance peaks. Furthermore, 
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if solar abundances Az Q are assumed, the equivalent hydrogen column density can be calcu- 
lated separately from each ion using the relation: 



N- 

1 y ion 



' H 



fion (Jimax ) Az^ 



(8) 



and placed at the position of £ max on an A#(log£) plot. For this we employed the XSTAR 



code (IKallman fc Krolikl|l995l ) version 2.1kn3 to calculate /j on (log£) using the continuum de- 



rived in extrapolated to t he range of 1 - 1 000 R vdberg. The results for IRAS 13349+2438 



using solar abundances from lAsplund et al.l (120051 ) are presented in Fig. |H The ions per 



taining to each element in the figure are connected by straight lines to guide the eye. The 
fact that the connecting lines of Fe, Si, Mg, O and Ne largely overlap on this plot implies 
that the relative abundances of these elements do not deviate much from the solar values. 
On the other hand, the relative abundance of N and Na appears to be significantly above 
solar, although there is only one N and one Na ion to show this result. The errors in Fig. H] 
are those propagated from the N ion symmetrical uncertainties as explained in the previous 
section. 

The distribution presented in Fig. H]is only a first approximation for the actual AMD. 
For the real AMD, we need to find a distribution dN H /d(\og^) that after integration (eq. [7j) 
will produce all of the measured ionic column densities (Tabled]). In fitting for the AMD one 
must take into account the full dependence of fi on on £. We assume all charge states see the 
same ionizing spectrum. This is justified by the absence of significant bound-free absorption 
edges in the spectrum. All elements are expected to reflect the same AMD distribution, due 
to the assumption that they all reside in the same gas. Iron however, has a special role as it 
covers almost five orders of magnitude in £, more than any other element. Nonetheless, in 
order to improve the AMD fit, we can use other elements as well. In order to incorporate the 
non-Fe ionic column densities without assuming anything about the elemental abundances, 
for non-Fe ions we use column density ratios. Using a reference ion for each element (the 
one with the smallest symmetrical errors), these ratios can be written as: 



N ion M43|)/Uiog£MiogO 



N ion-reference A z J ^(lo^) f ion-reference 



(9) 



Since there is no a-priori physical argument for a specific functional form for the AMD, 
we choose to use a simple staircase function (Fig. [5]). The main advantage of this form 
is the well-localized easily-calculated errors. A similar approach was used for emission by 



Nordon et al.l (120061 ) 
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The fit is performed with our Monte-Carlo evolution program in which the only free 
parameters are the values of the AMD in each log£ bin. The program receives starting values 
for each bin, then generates 5000 other AMD solutions randomly in the multi-dimensional 
region around the starting point. For each solution, x 2 is calculated and the lowest-^ 2 
solution is subsequently considered as the new starting point. This process is repeated until 
X 2 can not improve any more. In our tests, the program always converged to the same 
best- AMD solution for arbitrary starting values within five orders of magnitude (higher or 
lower) of the final result. This gives us confidence that the fit finds the global minimum. As 
this is an ill-defined inversion problem, the AMD binning is driven by the correlated errors. 
At first, we use 10 bins in the range of log£ = -1.25 to 3.75 (c.g.s.), assuming the AMD is 
zero elsewhere. We then combine bins so that the range and number of bins are optimized so 
that AMD bin values are inconsistent with zero (to within la), except in broad regions of £ 
in which the AMD is persistently zero. The final number of bins (i.e., the AMD resolution) 
depends on the specific absorber and on the quality of the data, which can limit the bin size. 

The AMD fit is obtained by the x 2 minimalization technique with respect to all of the 
measured Fe N ion values (eq. Ej) and N ion ratios of at least one other element (eq. [9]). In 
the calculation of x 2 , the real asymmetrical errors are used so that when the best-fit AMD 
overestimates an ionic column density, the upper limit error is used and vice versa. For non- 
Fe ions, the same rule is applied to N ion , while symmetrical errors are used for N ion _ re f erence 
(eq.EJ). 

The AMD errors are calculated in each bin including correlations between bins. The 
AMD in a given bin is varied from its best-fit value and the whole distribution is refitted. 
This procedure is repeated until Ax 2 = 1. The fact that changes in the AMD in one bin 
can be compensated by varying the AMD in other bins dominates the AMD uncertainties. 
This is what limits the number of bins and the AMD resolution. AMD in neighboring, 
excessively narrow bins can not be distinguished by the data, i.e., different distributions 
produce the measured Ni on values to within the errors. A meaningful quantity is the integral 
of the AMD up to £ (as a function of £) since during integration correlated errors cancel out, 
which explains the small errors on the integrated column density (Fig. [51 bottom panel). 

The next step is to estimate the element abundances relative to Fe. This is done by 
assuming that the other elements are distributed in log£ similar to Fe, or in other words, 
that the shape of the AMD (e.g., Fig. [5]) is manifested in the column density distribution of 
each and every element. This is what one would expect for a chemically uniform absorber. 
Indeed, the overlap of elements in Fig. @] seems to support this assumption. In order to 
obtain the abundances, we used the AMD distribution, by varying the AMD normalization 
to best-fit the column densities of each ion. This single-parameter fit yields for each ion 
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the relative abundance Az/Ap e . The errors on these abundances can be calculated from 
the Ax 2 = 1 requirement, again using the asymmetrical column density errors, but also the 
uncertainties on the AMD distribution itself. We then average this result for all ions from 
the same element to obtain our best abundance estimate. The errors on the abundances are 
calculated by averaging the upward and downward confidence limits separately. 



5.2. Results 



The best-fit AMD for the absorber in IRAS 13349+2438 is presented in Fig. El The 
integrated column density for the absorber of IRAS 13349+2438 is also presented in Fig. 
in the lower panel. This AMD was obtained using 21 ions of Fe and 9 ions of Si. The 
integrated AMD of the absorber in IRAS 13349+2438 (Fig. Qp gives a t o tal co lumn density 
of (1.2+0.3) xlO 22 cm" 2 compared with (1-4) xlO 22 cm" 2 in lsako et all (bpOlh . The AMD 



features a statistically significant minimum at 0.75 < log£ < 1. 75 (c.g.s. units), which 
corresponds to 4.5 < logT < 5 (K) and which was noted also by ISako et al.l (120011 ). This 
minimum is a manifestation of the relatively low ionic column densities observed for the ions 
Fe +n - Fe +15 seen in Fig. H and in Table HJ 

In order to check that the observed minimum in the AMD of IRAS 13349+2438 is not a 
result of the S/N, we want to apply the same method to another ionized absorber with a wide 
range of ionization. The natura l choic e is the 900 ks HETGS spectrum of NGC 3783. We use 
the spectrum from lKaspi et al.l (120021 . Fig. 5), and the same method described above. The 
results are presented in Figs. [6] and [7] in the same format as for IRAS 13349+2438 (Figs. H] 
and El). The AMD of NGC 3783 was obtained using 23 ions of Fe and 4 ions of Ne. Already 
in Fig. El low column densities for the Fe +12 - Fe +16 and Si +9 - Si +n charge states can be 
seen around log£ ~ 1-2 (c.g.s. units). Indeed, in the staircase AMD reconstruction for 
NGC 3783 presented in Fig. H, we find a minimum at 0.75 < log£ < 1.75 (c.g.s. units), 
which corresponds to 4.5 < logT < 5 (K). This minimum is surprisingly similar to that of 
IRAS 13349+2438, but is more statistically significant and thus better constra ined in £. We 



note that a similar minimum has been recently observed in NGC 7469 as well (IBlustin et al. 



20071 ). 



In Figs. [6] and [71 we include comparis ons with the two-com ponent model of lKrongold et al 



(120031 ) and the three-component model of lNetzer et all (120031 ) for the same NGC 3783 HETGS 
data set. The ionization components of these papers are placed on the plot according to T 
rather than £ as different models produce different relations between temperature and £. The 
bin size, is set arbitrarily to Alog£ = 0.5 (c.g.s. units). The column density of each compo- 
nent in those papers then determines the area (height) of its AMD contribution in Fig. [7J 
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It can be seen that all models avoid the region of the AMD minimum between log£ ~ 1 
and 2 (c.g.s. units). On the other hand, only the present method, that allows for a complete 
distribution, accounts for the full range of absorbing gas from log£ = —1.25 to log£ = 3.75 
(c.g.s. units). We obtain a total (integrated) N H value of (4.1±0. 7)xl0 22 cm" 2 comp ared 



wi th (3.8±1.0)xl0 22 cm; 2 and (2.0±0.8) xlO 22 cm" 2 obtained bv iNetzer et all (120031 ) and 



by lKrongold et al.l (120031 ). respectively. Note, that we use slightly different solar abundances 



than the other authors do. 

Finally, the results for the chemical abundances obtained for both IRAS 13349+2438 
and NGC 3783 using the method described in the previous section are listed in Table [2j We 
include in the table only elements with at least two significant ion detections. In both AGNs 
the Ne/Fe, Mg/Fe and O/Fe abundance ratios are slightly above or at about their solar 
value. The Si/Fe ratio is above solar. The S/Fe ratio can be measured only for NGC 3783 
where it is at its solar value. In general we can say that in both AGNs the metal abundances 
are approximately consistent with solar values. 



5.3. Discussion: Observing Thermal Instability 

The extent to which we can describe the low ionization region of the AMD strongly 
depends on the photoionization balance, which is highly uncertain at low-£ for two main 
reasons. First, the uncertainties in the dielectronic recombination rates of the low charge 
states of Fe at low (photoionized) te mperat u res ar e well known t o be seve r ely un derestimated 



in all codes including XSTAR; see iNetzerl (120041 ) and recently iBadnelll (120061 ). Correcting 



these would shift the ionization balance at a given £ value towards lower charge states, 
which would result in a more compact AMD not extending as far on the low-£ end. The 
other effect that could drastically change our low-£ AMD is photoionization of the low charge 
states by the EUV and UV continuum, which is only poorly constrained. For example, a 
softer continuum than we assumed here (see §3.ip . would allow for more highly ionized M- 
shell charge states at a given £ value. Consequently, it is important to understand that the 
exact position of the low-f peak of the AMD for both IRAS 13349+2438 and NGC 3783 is 
only poorly constrained, and will depend on the combined correction for the DR rates and 
for the soft continuum. Conversely, the total column density in this component as well as 
the existence of a minimum in the AMD are robust. 

One explanation that comes to mind for the minimum observed in the AMD of both 
IRAS 13349+2438 and NGC 3783 is thermal instability. If photoionized gas is thermally 
unstable at these temperatures, it would provide a natural explanation for the absence of 
absorption and for the low column densities of the charge states that form primarily at these 
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temperatures. Note that our AMD binning is driven by the requirement for meaningful 
errors. Consequently, we can not determine where the alleged instability occurs to better 
accuracy than the £ and T regimes given above, as narrower bins will inevitably result in 
excessive AMD uncertainties. The resulting AMD distributions presented in Figs. and [7] 
along with their similar minima are reminiscent of what is widely referred to in the l iterature 
as the two-phase models of gas in photoionization equilibrium (IKrolik et al.lll98ll ). These 
models have been invoked for AGN most importantly to allow for the confinement of low- 
T high-density clouds by high-T diffused gas. It is important to stress that the currently 
derived AMD is not a model. It represents a direct measurement of the ionization structure 
of the absorber that is more complete and more detailed than a multi iso-£ fit. To that 
end, the observed AMD minima provide direct evidence for a thermal unstable temperature 
region that can not be obtained from iso-£ fits. 

The observed minimum or two-phase structure can also be ascribed to two geometrically 
distinct regions along the line of sight, a high ionization region and a low ionization region, 
both which have their own gaussian like AMD distribution. However, there are good reasons 
to believe that gas at 4.5 < logT < 5 (K) would be unstable as th e cooling function A(T) is 
ge nerally decreasing in this temperature regime (IKrolik et al.lll98ll ). Additionally, the model 



of Netzer et a. 



The model of 



(2003 , Fig. 12) 



Goncalves et al 



or N GC 3783 shows that this region is marginally unstable. 



(120061 ) for NGC 3783 avoids temperature zones, which could 
be a result of thermal instability. Interestingly, the avoided region at 4.5 < logT < 5 (K) 
is consistent with the present measurement. However, there is another region of avoidance 
in that model at 5 .3 < logT < 6 (K) that our measurement does not show. The model 



of IKrongold et al. 



both Netzer et al. 



(2003) has a very wide uns t able region at 5.3 < logT < 7.5. However, 
(120031 ) and IKrongold et al.l (120031 ) find an absorption component in this 
temperature regime, as we do. As the regions of instability depend strongly on unobserved 
regions of the ionizing spectra, the theoretical stability analysis models have to be treated 
with caution. 



6. CONCLUSIONS 

We have analyzed the thermal and chemical structure of the ionized outflow in the quasar 
IRAS 13349+2438. Using a new AMD reconstruction method, we measure the distribution 
of column density in the outflow as a function of £. We find a double-peaked distribution 
with a significant minimum at 0.75 < log£ < 1.75 (c.g.s. units), which corresponds to 
temperatures of 4.5 < logT < 5 (K). Using a comparison method of the AMD derived for 
different elements we are able to estimate the relative chemical abundances in the outflow. 
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We find that all of the abundances relative to Fe are more or less solar. For comparison, we 
applied the same analysis to the excellent-quality X-ray spectrum of NGC 3783. We find a 
similar minimum in the AMD and rather similar abundances suggesting perhaps that these 
could be common features in AGN outflows. We believe this minimum is due to thermal 
instability. 

This research was supported by grant #28/03 from the Israel Science Foundation, by a 
grant from the Asher Space Research Institute at the Technion. S.K. acknowledges support 
of the Zeff fellowship at the Technion. 
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Table 1. Best-fit column densities for ions detected in the 2004 HETGS spectrum of 
IRAS 13349+2438 and comparison to the 2000 RGS observation jSako et alJboOlh . 



Ion HETGS RGS Ion HETGS RGS 

Column Density Column Density Column Density Column Density 
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Table 2. Relative abundances with respect to Fe in the ioniz ed X-ray absorbers o f 
IRAS 13349+2438 and NGC 3783, compared with solar ratios JAsplund et al.ll2005| ) 



Element 


Solar Ratios 
(A z /A Fe ) Q 


IRAS 13349+2438 

(A z /A Fe ) 


NGC 3783 

(A z /A Fe ) 
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{k Z /k Fe )/{k Z /k Fe ) Q 


NGC 3783 

(A z /A Fe )/(A z /A Fe ) 
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5 +ai 
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a Number in parenthesis indicate the number of ions used for each element, respectively for IRAS 13349+2438 
and for NGC 3783 in the determination of the abundances. 
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Fig. 1 — The present HETGS spectrum of IRAS 13349+2438 (in green) observed on 22-26 
February, 2004 compared to the RGS spectrum (in black) observed on 19-20 June, 2000. 
Despite the ~ 4 years difference, the flux level of the two observations appears almost 
unchanged, especially the hard X-Ray region and up to 17 A. There is noticeable variability 
beyond 17 A. Two O line locations are marked to demonstrate the absorption variability of 
oxygen. 
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Fig. 2. — Profile of the strong, isolated Fe +16 absorption line at 15.01 A plotted in velocity 
space with the best-fit profile in red. The global best-fit position of —300 km s _1 (blue-shift) 
and best-fit broadening corresponding to FWHM of 1500 km s -1 (tWb = 640 km s _1 ) are 
demonstrated. 
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Fig. 3. — Chandra HETGS spectrum of IRAS 13349+2438 corrected for cosmological redshift 
(z = 0.10764). The red line is the best-fit model in which a blue-shift velocity of -300 km s _1 
was applied to all ions. Ions responsible for the strongest lines and blends are marked above 
the data. 
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Fig. 4.— Equivalent N H distribution (eq. ED obtained fo r the IRAS 13349+24 38 outflow 
assuming ions form at £, max and assuming solar abundances (jAsplund et al.l 120051 1 . Lines are 
drawn between data points just to guide the eye. Vertical offsets between elements indicate 
deviations from solar abundances. The corresponding temperature scale, obtained from the 
XSTAR computation is shown at the top of the figure. 
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Fig. 5.— AMD of IRAS 13349+2438 scaled by the solar Fe/H abundance 2.82 xlO" 5 . The 
corresponding temperature scale, obtained from the XSTAR computation is shown at the 
top of the figure. The third bin value is zero, so only an upper limit error was calculated 
Note the minimum at 0.75 < log£ < 1.75 (c.g.s. units) corresponding to 4.5 < logT < 5 
(K). The integrated column density up to £ is shown in the lower panel yielding a total of 
N H = (1.2 ± 0.3) x 10 22 cm" 2 . 
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Fig. 6. — Equivalent N H distribution (eq.[8]) o btained for the NGC 3783 outflow assuming ions 
form at £ max and assuming solar abundances (jAsplund et al.ll2005l ) . Lines are drawn between 
data points just to guide the eye. Vertical offsets between elements indicate deviations 
from solar abundances. The corresponding tempera ture scale, obtained from the XSTAR 
computation is shown at the top of the figu re. The iNetzer et al.1 (120031 ) three component 
model results and the iKrongold et al.l (120031 ) two component model results are plotted for 
comparison. 
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Fig. 7.— AMD of NGC 3783 scaled by the solar Fe/H abundance of 2.8 x 1(T 5 . The 
corresponding temperature scale, obtained from the XSTAR computation is shown at the 
top of the figure. Note the minimum at 0.75 < log£ < 1.75 (c.g.s. units) corresponding 
to 4.5 < logT < 5 (K), a behavior similar to that of IRAS 13349+2438, but better con- 
strained here. The fourth bin value is zero so only an upper limit error is calculated. 
The integrated column density up to £ is shown in the lower panel yielding a total of 
Nh = (4.1 ± 0.7) x 10 22 cm -2 . Green dashed and cyan do t ted li nes repr esent the three 



and t wo- components results, respectively, of iNetzer et al.l (120031 ) and of iKrongold et al. 



(120031 ) where each component is broadened arbitrarily to Alog£ = 0.5 (c.g.s. units) bins. 



